Only for GSM (c.f. Yoshimura, K. and M. Kanamitsu, Dynamical global downscaling of global reanalysis, submitted to Mon. Wea. Rev.)
abstract:
program nudges the Fourier coefficients of forecasted run ("fcst data")
with those of reanalysis (or other external) run (called "base data").
program starts with the conversion of spectrum base data to
fourier series. this conversion is taken place in every 6 hours
(nesthour). the base data is linearly interpolated for the time
in between. the spectra of fcst data are also converted to
fourier series.
two kinds of fourier series are nudged with a weighting constant as,
fnew=(1/(a+1))*(fold+a*fbase)
for those with longer wavelength than lnest (e.g. 1000km)
(for those with shorter, a=0).
a can be variable, too (attenuating function).
all of the variables, t,q,u,v,ps,gz can be nudged, but practically,
ps,u,v and t should be nudged. instead, zonal q (the first coefficient
of the fourier series) is replacable to that of the base data.
default setting:
u,v,t,ps : nudging
q : zonal average replace
n-scale : 2000 km
n-weight : attenuating function (indaa=1)
Specification
New codes/scripts/data
src/fcst/gloopn.F: main routine for scale-selective nudging.
>> There are some versions for different nudging methods. Namely,
gloopn_Div.F : Only divergence component will be nudged.
gloopn_Vor.F : Only vorticity component will be nudged.
gloopn_diur.F : Diurnal component will be eliminated
>> >>> Monthly averaged 6-hourly base-files, "badiur.YYYYMMHH" (00Z, 06Z, 12Z, 18Z) should be created beforehand)
src/fcst/basefnam.F: getting a file name "base.YYYYMMDDHH"
Modified codes/scripts/data
src/fcst/gsmstep.F: call gloopn (after "filtr2")
src/fcst/Makefile.in
def/misc.h
configure-model
install
Automatic Installation
add --enable-gsmdamp
./install --model=gsm --enable-gsmdamp
Step by step
1. obtain sources
cvs co src/fcst/gloopn.F src/fcst/basefnam.F src/fcst/gsmstep.F
2. make a binary.
cd gsm
./configure-model --enable-gsmdamp
cd src/fcst
# edit and change parameters (lnest, aa) if necessary.
make clean
make
3. Get re-analysis data (sigma-level data).
R2/R1 sigma files (T62) are in ECPC data server (tsuzumi.ucsd.edu).
sftp dataman@tsuzumi.ucsd.edu
# get appropriate files
4. Change resolution (e.g. T62->T248) and save as "base.YYYYMMDDHH"
Do scnvrt and chgr.
The base.YYYYMMDDHH files should be saved or linked in the executive directory.
General Information
abstract: program nudges the Fourier coefficients of forecasted run ("fcst data") with those of reanalysis (or other external) run (called "base data"). program starts with the conversion of spectrum base data to fourier series. this conversion is taken place in every 6 hours (nesthour). the base data is linearly interpolated for the time in between. the spectra of fcst data are also converted to fourier series. two kinds of fourier series are nudged with a weighting constant as, fnew=(1/(a+1))*(fold+a*fbase) for those with longer wavelength than lnest (e.g. 1000km) (for those with shorter, a=0). a can be variable, too (attenuating function). all of the variables, t,q,u,v,ps,gz can be nudged, but practically, ps,u,v and t should be nudged. instead, zonal q (the first coefficient of the fourier series) is replacable to that of the base data. default setting: u,v,t,ps : nudging q : zonal average replace n-scale : 2000 km n-weight : attenuating function (indaa=1)Specification
- src/fcst/gloopn.F: main routine for scale-selective nudging.
>> There are some versions for different nudging methods. Namely,- gloopn_Div.F : Only divergence component will be nudged.
- gloopn_Vor.F : Only vorticity component will be nudged.
- gloopn_diur.F : Diurnal component will be eliminated
>> >>> Monthly averaged 6-hourly base-files, "badiur.YYYYMMHH" (00Z, 06Z, 12Z, 18Z) should be created beforehand)Automatic Installation
Step by step
1. obtain sources2. make a binary.
3. Get re-analysis data (sigma-level data).
- R2/R1 sigma files (T62) are in ECPC data server (tsuzumi.ucsd.edu).
4. Change resolution (e.g. T62->T248) and save as "base.YYYYMMDDHH"- Do scnvrt and chgr.
- The base.YYYYMMDDHH files should be saved or linked in the executive directory.
in shell script, for example,export INPUT_RESOLUTION=t62k28 SCN=/home/work5/kei/GSM248/libs/etc/scnvrt CHGR=/home/work5/kei/GSM248/gruns/runscr/chgr INFILE=${INDIR}/sig.anl.${YEAR}${MONTH}${DAY}${TIME}.ieee TMPBASE=tmp.base OUTFILE=${OUTDIR}/base.${YEAR}${MONTH}${DAY}${TIME} TMPFILE=tmp.sig ${SCN} gsm osu1 62:192:94:28 -1:-1:-1:-1:-1 ieee:bin ${INFILE}:dummy ${TMPBASE}:dummy ${CHGR} ${TMPBASE} dummy ${OUTFILE} dummy osu15. modify the main script as necessary and execute.Notes